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Pinning and creation of vortices in superconducting films by a magnetic dipole 
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Interactions between vortices in planar superconducting films and a point magnetic dipole placed 
outside the film, and the creation of vortices by the dipole, are studied in the London limit. The 
exact solution of London equations for films of arbritrary thickness with a generic distribution of 
vortex lines, curved or straight, is obtained by generalizing the results reported by the author and 
E.H. Brandt ( Phys. Rev. B 61, 6370 (2000)) for films without the dipole. From this solution 
the total energy of the vortex-dipole system is obtained as a functional of the vortex distribution. 
The vortex configurations created by the dipole minimize the energy functional. It is shown that 
the vortex-dipole interaction energy is given by — m • b^^^, where m is the dipole strength and 
b^^^ is the magnetic field of the vortices at the dipole position, and that it can also be written in 
terms of a magnetic pinning potential acting on the vortices. The properties of this potential are 
studied in detail. Vortex configurations created by the dipole on films of thickness comparable to 
the penetration depth are obtained by discretizing the exact London theory results on a cubic lattice 
and minimizing the energy functional using a numerical algorithm based on simulated annealing. 
These configurations are found to consist, in general, of curved vortex lines and vortex loops. 

PACS numbers: 74.25.Ha, 74.25.Qt 



I. INTRODUCTION 

The experimental study of artificial superconductor- 
ferromagnet systems has received a great deal of atten- 
tion lately. The magnetic, superconducting, and trans- 
port properties of a great variety of such systems, mostly 
superconducting films with arrays of magnetic dots or 
anti-dots placed in its vicinity, have been reported in the 
literature 0, S B • The main interest in this type of 
system is to enhance and modify pinning of vortices, and 
thereby increase the critical current and stabilize new 
vortex phases. 

From the theoretical point of view, these systems are 
also interesting because they alow theoretical predictions 
to be tested in detail experimentally, since the ferromag- 
netic structures responsible for vortex pinning are well 
characterized, and can be changed in a controlable way 
over a wide range of parametrs. The goal of theoreti- 
cal work is to understand how the presence of the fer- 
romagnet changes the equilibrium and non-equilibrium 
behavior of vortices. 

The interaction between vortices and the ferromagnet 
results from the action of the inhomogeneous magnetic 
field created by the ferromagnet on the superconduc- 
tor. This interaction is expected not only to pin vortices 
placed in the film by an aplied field, but also to create 
vortices, and even to destroy superconductivity in some 
regions of the sample. The theoretical problem posed 
by these systems is rather complex, since the equilib- 
rium vortex states in the absence of an applied field are 
non-trivial. The first problem that needs to be solved 
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is to calculate the vortex-ferromagnet interaction and to 
obtain the equilibrium vortex state resulting from the 
competition between it and vortex- vortex interactions. 
This paper solves this problem for a simple model con- 
sisting of a point magnetic dipole placed outside a planar 
superconducting film of arbitrary thickness in the Lon- 
don limit. First, the exact solutions of London equations 
for a film with a given distribution of vortices, consist- 
ing of a generic arrangement of straight or curved vortex 
lines, is obtained. Second, the total energy of the vortex- 
magnetic dipole system is calculated as a functional of 
the vortex distribution. The equilibrium vortex configu- 
rations generated by the magnetic dipole can then be ob- 
tained from the energy functional by minimizing it with 
respect to the vortex distribution. 

The main new results reported in this paper are: i) The 
proof that the vortex-magnetic dipole interaction energy 
is — m • b^^^, where m is the magnetic moment and b^^^ 
the magnetic field caused by the vortices at the dipole 
position, and that this energy can also be expressed in 
terms of a magnetic pinning potential for vortex lines of 
any shape, ii) The exact expression for the energy of the 
vortex-dipole system as a functional of the vortex distri- 
bution, iii) Approximate vortex configurations generated 
by the dipole in films of finite thickness. 

Earlier work on the above described model is restricted 
to the calculation of the interaction between a magnetic 
dipole and straight vortex lines. For a semi-infinite su- 
perconductor this interaction was obtained by Coffey |^ 
as — m-b^^^/2. Thin films were considered by Wei et. al. 

by Sasik and Hwa 7], and more recently by Erdin et 
al. 8]. The authors of Refs.|i,0 find the that the inter- 
action energy is also — m • b^^^, whereas those of Ref. |^ 
obtain the interaction as — m-b^^^/2 plus an extra term. 
In Refs. 0, II, m the creation of vortices by the dipole was 
investigated by minimizing the energy for simple config- 
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ur at ions. Recently, Milosevic, Yampolsky and Peeters 
10] obtained the energy of interaction between straight 
vortex lines and a point magnetic dipole for films af ar- 
bitrary thicknesses. These authors find that the interac- 
tion energy consists of two terms: one is — and 
the other is an interaction between the screening current 
generated by the dipole and the vortex. These authors 
also study the creation of vortices by the dipole by ex- 
amining the interaction energy of several configurations 
of straight vortex lines. 

This paper goes beyond these earlier results by con- 
sidering the interaction between the magnetic dipole and 
vortex lines of any shape. This is necessary in films that 
are not too thin, because the vortex-magnetic dipole in- 
teraction is limited to a distance of the order of the pen- 
etration depth from the film surface nearer to the dipole, 
whereas the vortex line energy grows with the film thick- 
ness. Thus, creation of straight vortex lines in thick films 
is energetically disfavored. 

To solve London equations for the vortex-dipole sys- 
tem this paper starts from the results obtained by the 
author and E.H. Brandt for films of arbitrary thick- 
ness without the dipole. In Ref. O the magnetic field 
and energy of an arbitrary vortex distribution in the film 
are obtained by solving London equations by the method 
of images. Here these results are generalized to include 
the magnetic dipole. Since London equations are linear, 
the total field of the vortex-dipole system is just the sum 
of the vortex field obtained in Ref. [ij with the field cre- 
ated by the magnetic dipole and the screening currents 
generated by it. From this result the total energy of the 
vortex-dipole system is obtained as the sum of the vortex 
energy in the absence of the dipole, and the vortex-dipole 
interaction energy — m • b^^^. The former, obtained in 
Ref. [Ill, is a quadratic functional of the vortex distri- 
bution, and is written here as the energy of interaction 
between the vortices in the film. The vortex-dipole inter- 
action energy is a linear functional of the vortex distri- 
bution, since in London theory the field b^^^ is a linear 
in the field sources. The functional coeficient of linearity 
is interpreted as the magnetic pinning potential. The ex- 
act expression for this potential is obtained here, and its 
dependence on the spatial coordinates and on the model 
parameters - magnetic moment strenght, position and 
orientaion, film thickness, and temperature - is studied 
in detail. It is also shown here that the vortex-dipole 
interaction energy is closely related to the screening cur- 
rent induced by the dipole: the change — m • Jb^^^ due 
to small deformation in the shape of the vortex lines is 
equal to the negative of the work done by Lorentz force 
of the screening current during the deformation. 

In the case of straight vortex lines the interaction en- 
ergy — m • b^^^, is found to be identical to that obtained 
by Milosevic et al [l^- Thus, with the exception of Ref. 
I3, the earlier results mentioned above for the energy of 
interaction between straight vortex lines and the mag- 
netic dipole agree one another and with the one obtained 
here. 



Minimization of the vortex-dipole system energy func- 
tional is not feasible in general because it involves infinite 
many degrees of freedom which are required to describe 
arbitrary configurations of curved vortex lines. Here the 
minimization is carried out approximately using the fol- 
lowing method. First, the exact London theory results 
are used to formulate a description of the vortex-dipole 
system on a cubic lattice. This description preserves the 
physics of London theory, and has the advantage that 
arbitrary configurations of vortex lines can be described 
by a finite number of variables. Second, the vortex-dipole 
system energy functional is written in terms of these vari- 
ables and minimized numerically, using simulated anneal- 
ing techniques. 

This paper is organized as follows. In Sec. ^ the ex- 
act solutions of London equations for the vortex-dipole 
system are obtained, and the total energy is calculated. 
In view of the diversity of formulas for the vortex-dipole 
interaction energy obtained by the earlier workers cited 
above, the energy calculation in Sec. ^ is carried out in 
detail. Two models for the magnetic dipole are consid- 
ered: a small current loop and a permanent dipole. In 
order to obtain the vortex-dipole interaction energy as 
— m-b^^*^ it is fundamental to use particular properties of 
the solutions of London equations reported in Ref. 11. In 
order to make the contact with Ref.jll easier, the present 
paper uses the same notation. The mathematical details 
of the calculations described in Sec. ^ are given in Ap- 
pendix^ and the relationship between the vortex-dipole 
interaction energy and the screening current is demon- 
strated in Appendix ini In Sec.|lll|the results of Ref.[ll| 
are used to write the total energy of the vortex-magnetic 
dipole system as a functional of the vortex distribution. 
First, in Sec. IIII Al the vortex energy in the absence of 
the dipole is written in terms of vortex-vortex interac- 
tions. Then, in Sec. IIII Bl the vortex-magentic dipole in- 
teraction energy is written in terms of a magnetic pinning 
potential and the dependence of this potential on the spa- 
tial coordinates and parameters of the model are studied 
in detail. Applications of these results of are considered 
in Sec. II VI First, the energy functional for straight vor- 
tex lines is obtained and compared with earlier results in 
Sec. IIV Al Then, minimization of the vortex-dipole sys- 
tem energy functional is discussed in Sec. IIVBI Finally, 
the conclusions of this paper are stated in Sec. 



II. VORTEX-MAGNETIC DIPOLE 
INTERACTION 

The film is assumed to be planar, with surfaces parallel 
to each other and to the x — y plane, and of thickness d 
occupying the region —d < z < {). The superconductor 
is isotropic, characterized by the penetration depth A. 
The magnetic dipole, m, is placed above the film at tq = 
(xo,^o,^o) = (ro±,^o) (^0 > 0). (Fig.IH). 

The magnetic field of the combined vortex-dipole sys- 
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FIG. 2: Physical interpretation of the vectorial vorticity dis- 
tribution. The normal n is parallel to v{y). 



FIG. 1: Superconducting film and magnetic dipole, m, at 
(ro± = 0,2;o). An example of vortex lines (full lines) and 
respective images (dashed lines). 



tern is written as 

b(-) = 

^j(out) _ 



^vac_ 



-b^ {-d<z<Q) 
b^^' {z <-d,z>Q), 



(1) 



where b^^"^ and b^^^ are, respectively, the vortex mag- 
netic fields inside the film and in vacuum. The fields bj^ 
and b^* are the dipole magnetic fields inside the film 
and in vacuum, respectively. The vortex fields where ob- 
tained in Ref. 11 by the method of images. According to 
it, any configuration of vortex lines in the film is char- 
acterized by a vectorial vorticity distribution u(y) or its 
Fourier transform i>'(k) defined as 



r e 



ikr 



(2) 



For vortex lines with vanishing core diameter, ^'(k) is 
given by a sum of line integrals along the vortices. 



z.(k) = Y^j 



d\j e' 



ikr. 



(3) 



The physical meaning of the vectorial vorticity distribu- 
tion is that the flux of i/(r) through an area perpendicular 
to it is an integer whose absolute value is the number of 
fiux quanta carried by the vortex line and the sign is that 
of the magnetic flux through the surface, as illustrated 
in Fig.ia 

Inside the film, the magnetic field b^^"^(r) satisfies the 
inhomogeneous London equation 



(4) 



Outside the film, assuming vacuum, the magnetic field, 
b^^*^ can be derived from a scalar potential that satisfies 
the Laplace equation 



bvac ^ _y^^ V^^ = 0. 



(5) 



The boundary conditions at the surfaces between the su- 
perconductor and the vacuum (2: = and z = —d) are 
that the perpendicular component of the current van- 
ishes and that the magnetic field is continuous. The 



method of images defines a vortex distribution in all 
space (—00 < z < 00 ), i>'^^(r), such that the current 
generated by it satisfies the boundary conditions, and 
that it coincides with the prescribed vorticity inside the 
film. In Ref. |il| it is shown that ly^^ consists of the vor- 
tex distribution u and its specular images at the film 
surfaces. This gives a periodic vortex distribution in 
the 2:-direction with period 2d. For the basic interval 
—d < z < d , L/^^(r) is given by 



i/^^(r) = 

u^l{x,y,z) = 
vl\x,y,z) = 



^(r), 

-iy±{x,y,-z), 
u^{x,y,-z), 



-d<z<0, 
0< z<d, 
0<z<d, (6) 



where ± stands for the vector component parallel to the 
x — y plane. An example is shown in Fig.^ The magnetic 
field inside the film is then 



j^stray 



(7) 



where b^^ is the field produced by the vortex distribution 
and its images, and is obtained by solving London equa- 
tion, Eq. in all space with i/^^{r) as the field source. 
The stray field inside the film, b^*^^^, is a solution of the 
homogeneous London equation. 

The magnetic field of the dipole inside the film, bj^, is 
a solution of the homogeneous London equation. Outside 
the film, b^* is the sum of the magnetic dipole field in 
the absence of the superconductor and the field of the 
screening current induced by the dipole. The boundary 
conditions are the same as those for the vortex fields. 
The magnetic field of the dipole is discussed in detail in 
Ref. [13 

The total energy of the vortex-dipole system, Et^ de- 
fined as the sum of the kinetic energy of the supercurrent 
in the film with magnetic energy inside and outside the 
film, can be written as 



Et = E]r 



(8) 



where Ei^ is the London energy of the supercurrent and of 
the magnetic field inside the film, and £^out is the energy 
of the vacuum magnetic field 



(out) 



d'^r± 
Stt L 



dz 



V X b(^") |2 + I b^^") |2] , 

d 

dz I b^^^*) I' . (9) 
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Now the couphng between the magnetic dipole and the 
superconductor is considered for two particular models 
for the dipole: a small current loop, and a point dipole. 

In the case where the magnetic dipole is a small current 
loop, the change in the total energy resulting from a small 
change in the vortex distribution, 5Et^ must be equal to 
the work done by the external source attached to the loop 
to keep the current constant during the change, SEext, 
that is 



SEt - SEe. 



where 



6Eext = 







(10) 



(11) 



where }ext is the current density in the loop, and (5a^^^, 
defined by Jb^^^ = V x ^a^^^, is the change in the vector 
potential of the b^^^ field. Expanding ^a^^^ around tq 
as 

fe""^(r) = fe"""(ro) + (r - tq) • V(5a""^(ro) 

it is straightforward to show that 

SEext = ni-Sh^^^{ro) . (12) 

It turns out that this is the only vortex-dipole interac- 
tion term in Eq. (|Tnj) . Other possible contributions to the 
vortex-dipole interaction would result from cross terms in 
SEt containing products of the vortex and dipole fields. 
It is shown in Appendix ^ that these terms vanish, so 
that SEt is the same as in the absence of the magnetic 
dipole. According to these results = — m-5b^^^(ro) 

can be interpreted as the change in the energy of inter- 
action between the vortices and the dipole. 



En, 



(13) 



The same vortex-magnetic dipole interaction energy is 
obtained if the magnetic dipole is modeled by a perma- 
nent point dipole. This model explores the well known 
analogy between magneto statics and electrostatics in 
current free regions of space 13]. According to it, b^* 
can be derived from a scalar potential, b^* = — V<l>rn 
that satisfies the Poisson equation 



-47r m • V(5(r — tq) 



(14) 



In this case the problem is identical to that of an electric 
dipole in an external field. The vortex-magnetic dipole 
interaction energy comes from the crossed term in £^out 
with b^* -b^^^, as shown in Appendix IXI. This approach 
is the same as that used in Refs. 6, 7, 10. In this case 
the vortex-magnetic dipole interaction energy is the work 
done to bring the dipole from infinity to its final position. 

To summarize, the total energy of the vortex-magnetic 
dipole system can be written as 

i?T = i?.-m-b™(ro)-im-b;„(ro), (15) 



where Ey is the energy of the vortex distribution alone. 
The last term in Eq. (|T^ is the energy of the dipole alone 
in the presence of the superconductor, b'j^(ro) being the 
field of the dipole screening current at the dipole posi- 
tion. In this paper this term is a constant, since m is not 
allowed to change. From here on this term is dropped. 

The vortex-dipole interaction energy, Eq. (|T3j) . can be 
generalized to any distribution of permanent magnetic 
dipoles placed outside the film and described by the mag- 
netization M. The result is 



M 



J 



(frM{r) ■ b™<^(r). 



(16) 



This expression is in agreement with a general formula 
of classical electrodynamics that expresses the interaction 
energy of a magnet in an external field in terms of the 
field that would exist in the absence the magnet Q. 
Here this field is b^^^(r). 



III. VORTEX-MAGNETIC DIPOLE ENERGY 
FUNCTIONAL 

Here the results of Ref. are used to express the total 
energy of the vortex-dipole system as a functional of the 
vortex vectorial distribution. 



A. vortex- vortex interactions 

It follows from Eqs.(32)-(35) and (22)-(28) of Ref. [Ill 
that the vortex energy, Ey, is a quadratic functional of 
the vectorial vortex distribution, which can be written as 



Ey/{eoX) 



dz dz' 



[g^ ( I rx - rl I ',z,z') (r^ ,z)-u^ (r^ ,z')^ 
g,s(\ r± -rx \]z,z')v,(v^,z)v,(v'^,z')] , (17) 

where eo = (^o/(47rA))^ is the basic scale for en- 
ergy/length. The dimensionless functions and Qzs de- 
scribe, respectively, the interactions between the compo- 
nents of the vorticity perpendicular and parallel to the z- 
axis: comes from vortex- vortex and vortex-image in- 
teractions, whereas Qzs has one contribution from vortex- 
vortex and vortex-image interactions, denoted Qz^ and 
another resulting from the energy of the stray and vac- 
uum fields, denoted Qsvi that is 

Qzs ; ^, z') = Qz{r±]z, z') + Qsv ; ^, z') (18) 

The functions Q^ and Qz are given by 



S^(±)(r±;^,^0 = ttA 



-t\z-z'\ 



e '^^coshrfz — , , coshrfz + + , , , 
- + (-) Z-T^, (19) 



sinh rd 



sinh rd 
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where r = \/k'^ -\- A~^, and the plus (minus) sign is for 
Qz ( Q±)' The first terms in the brackets in Eqs. (^HJ 
come from bulk vortex- vortex interactions, whereas the 
second and third terms come from the interactions be- 
tween the vortices and their images. These two functions 
need a short range cutoff to avoid unphysical divergencies 
at the vortex core . The interaction function resulting 
from the stray and vacuum fields is given by 



[coshr(z - 
[coshr(z - 

where 

/i 

c 



27r r (Pki_ e^^^-'^^ 
"T J (27r)2 k^TismhTdf^^^ ^ 
d) cosh t{z' + cZ) + cosh tz cosh tz'] + /2 x 
d) cosh tz' + cosh tz cosh t{z' + d)] } , (20) 



(fc_L + T)e^'* + {k±_ - T)e" 



2k^ 



{k^+Tfe' 



(21) 



The functions ^x, Qz and Qs 
to the transformations z < — 



are invariant with respect 
' z'^ and (z, z') {—z — 



(i, —z' — d). They depend only on d/\ and on ^/A, where 
^ is the vortex core radius. For thin films {d <C A) the 
term in is absent in Eq. (|17|) . and Qzs reduces to 
Pearl's result ^16j . 



B. magnetic pinning potential 

The vortex-dipole interaction energy is a linear func- 
tional of the vectorial vortex distribution, which can be 
written as 

/[^ dz 
d'^r± / — Uz{r±,z)Uym{^±,z) , 

(22) 

where Uyrn is the magnetic pinning potential. This result 
follows from Eq. and from the linear dependence 

of the scalar potential for the vacuum field, <l>, Eq. JSl, 
on the z-component of the vorticity obtained in Ref. [llj. 
Writing $ as 

^(r±,z) = Jd^r^J dz'v,{v\,z')x 

/C(|r^-rY |;z,/), (23) 

it follows from Eq. (|T^ that the magnetic pinning poten- 
tial is given by 

Uyrn{Y^,z) = • Vo/C(| ro± - rx \]Zq,z) (24) 
00 A 

The expression for the kernel /C follows from Eqs. (23), 
(25), (27), and (20) of Ref. Ill The result is 



/C(|r|;.o,.) = -(4.)^A^ (2.)^ 
(A:x + r)e^(^+^) - (A:^ - r)e-^^^+^) 




FIG. 3: Dependence of the magnetic pinning potential in the 
plane of the dipole on the coordinates x and z for a film with 
d = A. 



According to Eq. (|22|) (eoA)[/^^ is the energy of interac- 
tion of a vortex element parallel to the z-direction with 
the dipole. Note that the vortex-dipole interaction en- 
ergy, Eq. (|22|) . does not depend on the component of the 
vorticity perpendicular to the z-direction, u^. However, 
Vz and are not independent, since vortex lines form 
closed loops or lines that begin and terminate at the film 
surfaces {V ■ v = 0). The magnetic pinning potential, 
Uyrn^ depends only on the scaled variables d/A, zq/X and 
m/(j)oX. Since A depends on the temperature, Uym is 
temperature dependent, as shown next. 

Assuming that the dipole is located in the x — z plane 
and that to± = 0, the pinning potential can be written 



Uym{^±,z) 



rnx_ 
00 A 



cos (p ■ 



rriz dlC{r±; zp, z) 
00 A dzp 
S/C(rx; ^0, z) 



dr± 



(26) 



(25) 



where (p is the angle between and the x-axis. Simple 
analytical expressions Uym exist in two limiting situa- 
tions: large distances and thin films. The behavior of at 
large distances follows from the k± limit in Eq. (|25|) . 
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The result is 



-Stt 



{ 



cosh((z + d)/A) 
sinh (d/X) 

rrix r^X^ 
XT ^7^~ — 2Ts72 



} . (27) 



This expressions is vahd for \/~{T^^^r~ZQ) ^ A in films 
that are not too thin, and in thin films {d <C A) 
for ^{r\ + zl) > A = 2X^/d. In thin films for 
+ ^o) < is given by 



[/^^(rx,^;) -47r{{ 
cos 



1 



■[1 



^0 (ri 

^0 



ziyi^ 
]}. 



(28) 



Note that for thin films U^rn^ Eq. (|^ . is independent of 
the temperature, since A drops out. For films of finite 
thickness, and at short distances, [/^^ has to be calcu- 
lated numerically. Some results for [/^^ on the plane of 
the dipole are shown in Figs. [2l and ^ For m perpen- 
dicular to the film surfaces [rrix = 0), t/^rn is invariant 
under rotations around the 2;-axis, so that the graphs in 
Figs. El and left panels, represent the dependence of 
[/^rn on r^. For m parallel to the film surfaces (rriz = 0), 
the dependence of Uyrn on r± in a plane rotated by cp 
around the z-axis with respect to the x — z plane is just 
that shown in Figs. E^nd O right panels, multiplied by 
cosip (Eq. (|2S|)). These results show that the magnetic 
pinning potential penetrates a distance ~ A into the film, 
and that its range parallel to the film surfaces is a few A. 
The temperature dependence of Uym^ shown in Fig. [51 is 
that increasing the temperature increases both the range 
and the absolute value of Uyrn- 




E 
> 

Z) 


W 2^/^=0.25 







FIG. 4: Dependence of the magnetic pinning potential in the 
plane of the dipole for a film with d = A on the coordinate x 
at 2; = (full line) and z = —d (dashed line) for two values of 
the dipole height zq. Left panel: rux = 0, niz = 0oA. Right 
panel: rUx — 0oA, ruz — 0. 



IV. APPLICATIONS 

In this section the results derived above are applied to 
study the equilibrium states of vortices induced by the 
magnetic dipole. To obtain these states in the absence 
of an applied magnetic field it is necessary to minimize 
the energy Et, Eq. (|T5|) . with respect to the vorticity 
distribution i/. The equilibrium state, in the absence of 
an externally applied field, will contain vortices if there 
is a configuration of vortex lines for which Et = Ey -\- 
Eyrn is a minimum and negative. Minima with Et > 
are met ast able, since for the film without vortices Ey = 
Eym = 0- Minima with Et < can occur if there are 
configurations of vortex lines for which Eym is sufficiently 
negative to overcome the positive definite vortex-vortex 
interaction energy Ey. The shapes of the vortex lines 
generated by the dipole are dictated by the competition 
between Ey and Eym- The spatial inhomogeneity of the 
magnetic pinning potential also plays an important role 
in determining the shapes, as shown in Sec. IIVBI 

According to the discussion in Sec. lIIIl the equilibrium 
states depend only on the scaled parameters d/A, zq/X, 
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FIG. 5: Dependence of the magnetic pinning potential in the 
plane of the dipole, and at ^ = 0, on the coordinate x for a 
film with d = 2A at different temperatures, defined by A/A(0). 
Dipole height zo = A. Left panel: rUx = 0, rUz = (/)oA. Right 
panel: rrix — (/>oA, rUz — 0. 
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and m/^o-^^ and on the temperature. If thermal vortex 
fluctuations are neglected, which is appropriate for low- 
Tc superconductors, the temperature dependence comes 
only from A. For a given film and magnetic dipole 
zo and m are fixed, but the parameters d/A, zq/X and 
m/(j)oX change with temperature, leading to nontrivial 
changes in the equilibrium state, as discussed in Sec. 
ITVRl 

The problem of minimizing Et with respect to the 
vortex distribution is, in general, a formidable task. The 
simplest case is that of straight vortex lines perpendicu- 
lar to the film surfaces, where the functional dependence 
of Et on the vortex distribution can be described by a 
finite number of degrees of freedom, as discussed in Sec. 
IIV Al For films of finite thickness the vortex lines are 
not in general straight, and Et depends on infinite many 
degrees of freedom required to describe arbitrary config- 
urations of curved vortex lines. The minimization of Et 
can only be carried approximately, by reducing the de- 
grees of freedom to a discrete set. A method to do this 
is discussed in Sec. IIVBI 



A. straight vortex lines 

Here it is assumed here that the vortex lines in the 
film are straight and perpendicular to the film surfaces. 
In this case iy± = and Vz is given by 



is obtained from Eq. (|5SJ, and (|S3) as 



(29) 



(0 



Where qi = 0, ±1, ±2.. is the vorticity of the i-th vortex 
line and its position in the x — y plane. The energy 
is then 

Er/ieoX) = {E(ij) U,v{\ Ri - |) + 

QiUltHRi)} (30) 



where 



Uvv{\ ~ I) — 



and 



° dz 
A 



Uym (Ri , Z^ . 



(32) 



The interaction energy of a vortex line pair (R^ Rj)is 
(eoA) 2Uyy{\ R^— Rj |), and (eoA) Uyy{0) is the vortex hue 
self-energy. The vortex-vortex interaction energy Uyy is 
discussed in detail in Ref. O. 

The interaction energy of the vortex line with the mag- 
netic dipole is (eoA) UlZ^CRi). The expression for U^Z^, 



f/r (R) = -(4^) 

m 



(27r)2 



( 'k. 
(j)oX 0oA 



Ck±r 



(33) 



It is straightforward to shown that this expression is iden- 
tical to that obtained in Ref. 10.. It can be shown, us- 
ing the results of Appendix IbI that Ul^^{Il) can also be 
written as the negative of the work done by the Lorentz 
force of the dipole screening current to bring the vortex 
line from a position far from the dipole to R. In Ref. 
ilQi the vortex-dipole interaction energy is written as the 
sum of minus one-half of the Lorentz force work with 
-m • b^^^/2. 

The total energy Eq. (pUj) is a functional of the vor- 
ticities {qi) and positions (R^) of the vortex lines. In 
order to obtain the equilibrium vortex configurations Et 
has to be minimized with respect to these variables. The 
minimization involves only a finite number of variables, 
and can be carried out numerically with modest compu- 
tational resources. This minimization will be discussed 
elsewhere. 

The equilibrium vortex configurations created by the 
magnetic dipole will consist of straight vortex lines only 
for thin films (d <C A). For thick films d ^ X, the self en- 
ergy and the vortex- vortex interaction energy grow with 
the thickness of the film, (i, whereas the vortex-dipole in- 
teraction energy does not, since Uym is limited to region 
of depth A from the film surface closest to the dipole. 



B. lattice London model 

In this section an approximate method to obtain equi- 
librium vortex configurations induced by the magnetic 
dipole is presented. It is based on a discretization of the 
vortex degrees of freedom called lattice London model. 
This model was introduced several years ago to study 
vortex fluctuations in high-Tc superconductors 17]. In 
the present context the lattice London model is useful be- 
cause it requires only a finite number of degrees of free- 
dom to describe curved vortex line configurations, and 
because it preserves the essential physical ingredients of 
the vortex-dipole system in the London limit. 

The vortex distribution in the lattice London model 
is represented by integer variable placed on three- 
dimensional mesh with cubic unit cell of side a ~ ^, and 
subjected to periodic boundary conditions. At each lat- 
tice site there are three integers = 0, ±1, ±2..., one for 
each spatial direction ji = x^y^z. From the configuration 
of the variables at each lattice site, the configura- 
tion of vortex lines follow by associating arrows with the 



as shown in Fig. (SJa. Essentially, the lattice London 
model restricts the vorticity u to point in one of the x, 7/, z 
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FIG. 6: a) Graphical representation of integer vorticities at 
lattice point r. b )Vortex loops and vortex lines used in the 
numerical minimization method to generate vortex configura- 
tions. 

directions, and represents the flux of v through the 
face of the cubic unit cell perpendicular to the /i direction 
(Fig.Ela). For bulk and semi- infinite superconductors the 
lattice London model is an exact discretization of Lon- 
don theory on a cubic lattice, as shown in Ref. 0. For 
the problem under consideration here, the lattice London 
model is an approximation. It consists in replacing the 
film by a cubic mesh of lattice constant a ~ ^, where 
the vortices are defined as discussed above. The cubic 
mesh is subjected to periodic boundary conditions in the 
X and y directions. For the sake of simplicity, in what fol- 
lows it is assumed that the vortex lines generated by the 
dipole are in the plane of the dipole. This reduces the 
search for vortex line configurations to two-dimensions 
{x — z plane), so that Uy = 0. This assumption is justified 
for the equilibrium vortex configurations discussed here. 
The vortex-dipole system energy functional is obtained 
from the one derived above for the continuum London 
model as described next. The vortex- vortex interaction 
energy is taken as 

i?./M) = (f)2E,,, 

\^Tlx (^i 1 Zi^Tix {Xj , )r_|_ ( I Xj | , , Zj^ 

nz{xi, Zi)nz{xj , Zj)T zs{\ xi - Xj \\Zi, z^)} , (34) 
and the vortex-dipole interaction energy as, 

^t;m/(eoA) = nz{Xi,Zi)Tyrn{Xi,Zi) . (35) 

The functions Fx and V and are obtained from 
Qsv) and Uvrn-, respectively, by making the latter ones 
periodic in the lattice along the x and y directions. This 
consists in substituting = (k^^ ky) in Eqs. (|19|) . ((^n|) . 

and by K = [2a~^ sin {kxa/2)^ 2a~^ sin (/c^a/2)], and 




m =1.75 



FIG. 7: Equilibrium vortex configurations generated by a 
magnetic dipole perpendicular to the film surfaces, and lo- 
cated at zq\ from the point indicated by the arrow, d in 
units of A, rriz in units of 00 A. Mesh spacing a — 0.1 A. 



by replacing the integrals over in the same equations 
by sums over the reciprocal lattice of the cubic mesh. 
This procedure is justified by the exact results of Ref. 
17. 

To minimize the functional Et = Ey -\- Ey^^ Eqs. ((5^ 
and with respect to Ux and Uz simulated annealing 
is used , together with the following procedure to generate 
the vortex line configurations 17]. First, it is attempted 
to add vortex loops at every lattice site: square loops at 
sites not on the film surfaces and open loops at surface 
sites. Second, it attempted to add straight vortex lines, 
perpendicular to the film surfaces at every position Xi. 
This is illustrated in Fig. Elb. 

The numerical minimization of Et is carried out for 
a few parameter values with a/X = 0.1. The results are 
shown in Figs. [7| and El These figures can also be viewed 
as representing the evolution of the equilibrium vortex 
configurations for the film with d = 3.0A(0), at T = 0, 
with increasing temperature. The values of the scaled 
parameters d/A, zq/X and m/(j)oX are chosen so that the 
sequence of panels corresponds temperatures such that 
A/A(0) = 1, 3/2, 3,6. Note that in this case the mesh 
parameter in Figs. [7| and [SI cl ^ also changes with T, 
since ^ is temperature dependent. 

m perpendicular to the film surfaces: The equilibrium 
vortex configurations consist of vortex lines with a sin- 
gle flux quantum, that is | Ux | = | riz |= 1 along the 
lines, and with Uz of the same sign as (positive in 
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d=3.0 



m =7.5 



d=1.0 



m =2:5 



d=2.0 
z^=0.8 

m =5.0 



FIG. 8: Equilibrium vortex configurations generated by a 
magnetic dipole parallel to the film surfaces, and located at 
zqX from the point indicated by the arrow, d in units of A, 
TTix in units of 00 A. Mesh spacing a = 0.1 A. 



Fig.EJ. The vortex lines are curved, except for the thin- 
ner film with d = 0.5 A. The curvature results from the 
competition between vortex- vortex and vortex-dipole in- 
teractions. Closer to the z = surface the vortex-dipole 
interaction dominates, pulling the vortex lines towards 
the dipole, and keeping them perpendicular to the film 
surfaces. Deeper inside the film the vortex-dipole inter- 
action weakens and vortex- vortex repulsion separates the 
vortex lines apart. The creation of net vorticity by the 
magnetic dipole obtained here is particular to the iso- 
lated dipole. For a dipole array the net vorticity must 
vanish, even if the dipoles are far apart, due to the long 
range of vortex- vortex interaction in films 16]. 
m parallel to the film surfaces: The vortex configurations 
consist of half loops and pairs of vortex lines with oppo- 
site vorticity, both with a single flux quantum, as shown 
in Fig. [HI These configurations reflect the properties of 
the pinning potential for m parallel to the film surfaces 
shown in Figs. 01 and ^ For instance, in the case of 
the half loop for d = 3. OA, the resulting curve places the 
negative (positive) in regions of positive (negative) 
pinning. In the case of vortex lines, the nearly straight 
curve follows, essentially the pinning potential maxima 
for negative and the minima for positive n^. The re- 
sults shown in Fig. [51 also indicate how the vortices pen- 
etrate the film with increasing temperature, that is the 
transition from the Meissner state to the mixed state. 
The transition is continuous. Half loops penetrate and 



grow towards the interior of the film, eventually separat- 
ing in two lines of apposite vorticities. 

The equilibrium vortex configurations described above, 
for both orientations of m, are not expected to change if 
the restriction that they are in the plane of the dipole is 
lifted, since there would be no gain in the vortex-vortex 
or vortex-dipole interaction energies if some vortex lines 
were out of the x — z plane. 



V. CONCLUSION 

In conclusion then this paper solves exactly in the Lon- 
don limit the problem of vortices in a film of arbitrary 
thickness interacting with a point magnetic dipole out- 
side the film, and obtains from these solutions the vortex- 
dipole system energy as a functional of the vortex distri- 
bution, which is minimum for the equilibrium vortex con- 
figurations generated by the dipole. The results reported 
here can be generalized to any distribution of permanent 
magnetic moments placed outside the film. The numer- 
ical method to obtain equilibrium vortex configurations 
presented here is of general validity, and can be applied 
to three-dimensions and to other distributions of dipoles. 

The London limit used in this paper, is expected to 
break down for large values of magnetic dipole, because 
the inhomogeneous magnetic field created by it destroys 
superconductivity locally in the film. Roughly speak- 
ing, London theory is valid as long as the maximum field 
of the dipole at the film surface nearer to the it is less 
than the upper critical field, that is , tti/zq < 0o/(27r^^), 
or m/{(j)oX) < {zo/\Y{\/C)'^ /2tt. One indication of this 
break down is the appearance in the equilibrium vor- 
tex configurations of vortex lines separated by distances 
< 2^. In the equilibrium vortex configurations shown in 
Figs. [3 and[Hlthis occurs for m perpendicular to the film 
surfaces in the film with d = X, and for m parallel to the 
film surfaces in the films with d = X and d = 0.5A. In 
both cases the values of m are found to be in agreement 
with the condition stated above. Note, however that only 
in the immediate vicinity of the film surface the vortex 
lines separation is < 2^. Deeper inside the film the vor- 
tex lines are separated by distances larger than 2^. This 
can be interpreted as indicating that the regions where 
the vortex lines separation is < 2^ are normal, and that 
in the regions where the separation is larger than 2^ the 
vortex configurations obtained in the London limit are 
reasonable estimates. In the case of bulk superconduc- 
tors and films under applied magnetic fields, the inter- 
pretation along similar lines of London theory results for 
the vortex configurations generated by the field lead to a 
reasonable first order approximation to the vortex phase 
diagram. The same is believed to be true here. The 
London theory results described in this paper, and their 
generalization to distributions of dipoles, can be applied 
beyond their strict limits of validity to give a first order 
approximation to vortex behavior in these systems. 
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APPENDIX A: MATHEMATICAL DETAILS 







u(z) / 



Here some details of the derivations in Sec.^are given. 
First it is shown that all cross terms in Et vanish. 

When Eqs. (^) are substituted in Eq. there are 
several cross terms containing two distinct fields. Here 
it is shown that these terms vanish. The first step is to 
show that there are no cross terms in E^^ with b^^ and 
any homogeneous solution of London equation, denoted 
as b^. This term is 



Er = 



47r 



dz 



[A^Vxb^^-Vxbh + b^^-b^]. 



(Al) 



Using the identity 

V X b^^ V X b^ = 
V • [b^^ X (V X b^)] + b^i • [V X (V X b^)] , (A2) 

and the fact that b^ satisfies the homogeneous London 
equation, Ec can be written as 



Er = X^ 



l^z-[h-'x{Vxh^)]f_, . (A3) 



As shown in Ref. 11, b^^ at the film surfaces z = 0, —d 
points in the z-direction, so that the vector product in 
the integrand has no z-component, and Ec = 0. This 
argument eliminates the cross terms with b^^ and b®*^^^, 
and with b^^ and bj^. One cross term is left in Ei^ with 
the fields b^*^^^ and b^^^. It is shown next that in the case 
of a small current loop this term is canceled out by the 
cross term with b^^^ and b^* in £^out- Denoting these 
terms by Ei^c and £^outc, respectively, it follows that 



Eu 



J 4^ J-d 



dz X 



[A^V X b'^^'-'^y • V X bj^ + Yf^'""^ ■ bj^] , (A4) 



Eo^^t c. — 



47r 



/ dz ^ [ 

Jo J-( 



dz 



^vac _ j^out 



(A5) 



Using arguments similar to those leading to Eq. (| A3|) . 
Einc can be written as 

Einc = ^J dWi ■ [b^*-^"^ X (V X h^)] . (A6) 



FIG. 9: Full line: generic vortex line defined by u{z). Dashed 
line: vortex line after deformation Su{z). 



From London equation V x = aJ^/A^ (bj^ = V x aj^), 
so that 

^inc = ^ 1 d\^Z . K X b^*-^] 1° , . (A7) 

It is possible to write £^out c in a similar form, using V x 
= 0. After an integration by parts of Eq. (|A5|) the 
result is 



out -uvaci lO 



£^outc = -^ /dViZ-[ar xb 



(A8) 



It follows from the continuity of the fields and vector 
potentials at the film surfaces that Ei^c + ^outc = Cl- 
in the case of a permanent magnetic dipole there is 
only a partial cancelation, and Ei^c + ^outc = — m • 
b^^^(ro), as shown next. The energy £^outc is written in 
terms of scalar potentials using the identities 



= V<I>V<I>m 



-47r^m • V(5(r 



ro). 



Substituting in Eq. l|A5|l results 

£outc = -m-b--'=(ro)- 



(A9) 



The second term in Eq. (|A9|) cancels out Einc- This can 
be shown starting from Eq. (|A6|) . and using Eq.(22) of 
Ref. El' for h^^'^y, the homo geneous London equation for 
bj^, and continuity of the fields at the film surfaces. 



APPENDIX B: WORK DONE BY THE LORENTZ 
FORCE 

Here it is shown that the change in the vortex-dipole 
interaction energy SEym , Eq. ([T^ , equals the negative of 
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the work done on the vortex hnes by the Lorentz force of 
the screening current induced by the magnetic dipole. 

Consider the vortex hue running from one film surface 
to the other shown in Figs. El The equation describing 
this hne is 

t{z) = zz^u{z) . (Bl) 
The contribution of this hne to the vorticity is 



z/(k) = / dz[z^ 

-d d.z 



(B2) 



It is convenient here to work with Fourier transform in 
the X — y plane only. For the vorticity it is 



(B3) 



If the vortex line undergoes a small deformation, defined 
by 5\i{z)^ the change in the vorticity to first order is 



8v{\<i^,z) 



{[— zkx • 5\x{z)]z - 



d5\i{z) 

dz 



(B4) 



The corresponding change in the vortex-dipole interac- 
tion energy is 



dz 

(27r)2 J_, T 



(B5) 



When the vortex line is deformed by 5\i{z)^ the Lorentz 
force of the screening current induced by the magnetic 
dipole, j*^ = c/47rV x bj^, does the work 



[z+^!M].5u(^)e-*k^-"(-) . (B6) 

dz 



The screening current is perpendicular to the z-direction, 
and is given by (Ref. llQ) 



ckl ' 



Because both the screening current and 5\i{z) are parallel 
to the film surfaces, the term with du{z)/dz in Eqs. (|B6|) 
vanishes. Substituting Eq. (|B7|) in Eq. (|B6|) . and using the 
expression for Uym obtained in Sec. IIII Bl it follows that 
SEyrn = —SWl- This result can also be demonstrated for 
vortex lines that cannot be described by Eqs. (|B1|) . such 
as loops and lines with humps. 
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